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Abstract 



^ I The S = 1/2 XXZ spin chain with the staggered XY anisotropy 

^ • N N 

C , 11 = J^^{S^S^_^_i + S'^S^^i + AS^S^^i) — ^y^(— l)"'('5'^'S'j^+i — S^S^+i) 



is shown to possess gapless, Luttinger-hquid-hke phases in a wide range of its 
parameters: the XY-hke phase and spin nematic phases, the latter character- 
ized by a two-spin order parameter breaking translational and spin rotation 
symmetries. In the simplest, exactly solvable case A = 0, the spectrum re- 
mains gapless at arbitrary J and 6 and is described by two massless Majorana 
(real) fermions with different velocities v± = \ J ziz 6\. At \6\ < J the stag- 
gered XY anisotropy does not influence the ground state of the system (XY 
phase). At |5| > J, due to level crossing, a spin nematic state is realized, with 
TTii and tilt local symmetry of the xx and yy spin correlations. The spin 
correlation functions are calculated and the effect of thermally induced spin 
nematic ordering in the XY phase ("order from disorder") is discussed. 



The role of a finite A is studied in the hmiting cases \5\ <^ J and 

\S\ » J, using bosonization method. On the basis of the derived self-dual 

field-theoretical model, similar to the one recently proposed to treat two 

weakly coupled Luttinger chains, a confinement regime is described, in which 

fermionic excitations fail to be observable due to Luttinger-liquid effects. 

Kosterlitz-Thouless transitions to massive phases, driven by the anisotropy, 

are discussed. 
75.10.-b, 75.10.Jm 
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Typeset using REVT^ 



I. INTRODUCTION 



There has been considerable recent interest in the S = 1/2 quantum spin chains with 
alternating ferromagnetic and antiferromagnetic bondffl: 

^ = E[(S2,-S2,+l)-/3(S2,-l-S2,)]. (1) 

i 

The model (|1|) continuously connects a gapful phase of decoupled singlets {(3 = 0) with the 
limit of the 5* = 1 chain {j3 —>■ oo) and so far was mostly considered in the context of the 
Haldane gap problemH. 

In this paper, we would like to draw attention to another, in a sense, opposite aspect of 
the problem of bond-alternating spin chains. In the rotationally (SU(2)) invariant S = 1/2 
chain the formation of a massive phase appears as a result of bond alternation [jS ^ — 1)- 
On the other hand, as follows from the (Bethe ansatz) exact solution of the translationally 
invariant XYZ mo del, an Ising-like mass gap is always present in the excitation spectrum 
as long as continuous spin rotation symmetry remains fully broken. One then may ask: 
Are the requirements of translational and spin rotation symmetries robust for realization 
of a gapless, Luttinger-liquid low-energy behavior oi S = 1/2 spin chains? Can the bond 
alternation and full breakdown of spin rotational symmetry, each separately generating a 
gap, result in a gapless (critical) regime when acting together? As shown below (Sec. 2), 
among several unitarily equivalent possibilities, the XXZ spin chain with the staggered XY 
anisotropy 

H = Hxxz + Hs 

N N 
n n 

is a nontrivial example of such a system exhibiting gapless phases in a wide range of its 
parameters with rather interesting properties. 

In Sec. 3 we first consider the simplest, exactly solvable case of the XY bond-alternating 
chain (A = 0). The Jordan- Wigner transformation reduces it to a Id system of spinless 
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fermions with a half-filled band and a Cooper pairing with the total momentum equal vr. 
This model is also equivalent to two decoupled quantum Ising chains (i.e. Id Ising models 
in transverse magnetic fields) at criticality. Therefore, for any ratio 6/ J, the excitation 
spectrum is gapless, and elementary excitations are represented by two Majorana (real) 
fermions propagating with different velocities v± = |J± 6\a. At \6\ < J the staggered XY 
anisotropy does not influence the ground state of the system but reveals at the level of 
excitations, leading to splitting of the two gapless branches of the spectrum. This affects 
the time dependence of the spin correlation functions and thermodynamics. 

At \6\ > J, due to level crossing, the XY ground state is changed by another gapless 
phase - a spin nematic (SN) state, characterized by broken translational and spin rotational 
U(l) symmetry, but preserved time reversal invariance. In this phase, the local symmetry 
of the XX and yy spin correlations is changed by the TTii and tilt structures, respectively. 
This reflects the sign-alternating character of exchange couplings between neighboring spins 
at \6\ > J (see Eq.(0)). The new spin correlations (following the same power-law decay with 
distance as that in the XY phase) induce a two-spin SN ordering of the system described 
by the order parameter < vac\S!^S^^i — S^Sn+ilvac >~ (—1)" sgn 6. We also find a 
manifestation of Villain's "order from disorder"! - the effect of thermally induced SN ordering 
of the spins in the XY-like phase {\6\ < J). 

In Sec. 4 we extend consideration to the case of a finite A to study effects of a weak 
staggered XY anisotropy (|5| ^ J) in a Luttinger spin liquid. Using bosonization method, 
we show that the low-energy properties of the system (2) are described by a self-dual field- 
theoretical model, similar to the one recently proposed to describe two weakly coupled 
Luttinger chains!. In the range —\j\f2 < A < 1, the off-diagonal perturbation iY^, tending 
to split the velocities of fermionic excitations, is irrelevant due to the infrared catastrophe 
that eliminates single-fermion states from the low-energy part of the spectrum. This behavior 
is an example of Anderson's confinementi which has been extensively discussed recently 
in connection with the two-chain problem. In this regime, the system maintains basic 
properties of the gapless Luttinger liquid, slightly modified by thermally induced SN order. 
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At A < —\l\/2, a translationally invariant, two-axis anizotropy, generated in the second 
order in 5, becomes a relevant perturbation, and the system undergoes a Kosterhtz-Thouless 
transition to a Neel-ordered state with polarization axis x or y. 

By the 6 ^ J interchange symmetry, the above picture also holds in the opposite lim- 
iting case I (5 1 ^ J. Here one finds Luttinger-liquid-like SN phases and Kosterlitz-Thouless 
transitions to phases with ttii ordering of the spins, either along ct— or y— direction. 



II. CHOICE OF THE MODEL 

Consider a bond-alternating S = 1/2 quantum XYZ spin chain given by the Hamiltonian 

N 

^ = E E [Ja + i-irj:]s:s:^,, (3) 

n=l a=x,y,z 

The model can be mapped onto two coupled quantum Ising chains by a nonlocal unitary 
transformation, previously used by Kohmoto, den Nijs and Kadanofl^ in their study of the 
2D Ashkin- Teller model. This transformation is a special combination of spin rotation and 
duality transformation!. Its basic steps are: i) a global 7r/2 spin rotation, o"^ — ^ o"^, cr^ ^ 
-al, ii) duality transformation <_i/2<+i/2' «+i <+i/2> where {n+ 1/2} are 

the dual lattice sites; iii) relabeling the sites, j |(j + |)) introduction of two kinds 
of the Pauli matrices, a" and Tj!|.i/2) defined on the two sublattices of the dual lattice; iv) 
duality transformation of the r-spins only. As a result, spin operators 5*" of the staggered 
XYZ chain (2) are represented in terms of a and r matrices as products of order and disorder 
operators: 

OCX _ z nqx — rr^ t-2 

^^2i - '^i^i+1/2' ^'^2j+l - '^j+l^i+1/2' 

2S2j = ■> '^^Ij+l = '^|+l/2^/+l' 



where 



25"!^ — 'i<y'^<y^j^i/2'^'j'Tj+i/2i 2>S'|j_|_i — '^o"J+i/2^|+i'^/+i/2'^i+i' 



N/2 j 

^i+i/2= n ^r, ^/+i/2 = n^r- 

l=j+l 1=1 
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The transformed Hamiltonian 

H = Hfj + Hj- + H„T- 

N/2 

describes two coupled quantum Ising chains. The constants Aa = (1/4)(Jq, + J^), Ba = 
(l/4)(Ja-4). Following Ref. 10, it can be shown that this Hamiltonian determines the trans- 
fer matrix of the two-dimensional asymmetric Ashkin- Teller (AT) in the highly anisotropic, 
the so-called r-continuum limit. The asymmetric AT model represents two different 2D 
Ising models coupled by a four-spin interaction. 

To sort out those cases when the model (§) may show a critical (gapless) behavior, we 
shall treat the interchain ar coupling as a weak perturbation and then consider conditions 
when at least one of the two quantum Ising chains is critical. There are three such possibil- 
ities. 

1) ifo- is critical, while is not (or vice versa), implying that = ±By,Ay ^ i-B^- 
Using the Majorana fermion representation for each chain0, one finds one massless (i^o-) 
and one massive (-f^r) branch in the excitation spectrum. The interchain interaction H^r 
is irrelevant; its role is exhausted by renormalization of the mass and critical temperatute 
(or group velocity of the massless fermion). In this case, the critical behavior of the system 
belongs to the 2d Ising universality class. 

2) ifo- and Hr are identical and critical; e.g. = By = Ay = B^- If at the same 
time Az = Bz, the model (§) becomes equivalent to the exactly solved U(l)-symmetric, 
translationally invariant XXZ spin chainlll, related to the six-vertex mo delil. The inter- 
chain coupling H„j. is marginal: the system shows a Luttinger-liquid critical behaviour with 
continuously varying (coupling-dependent) critical exponents. 

3) H„ and Hr are both critical but different. This case occurs when 

A^ = ±By, Ay = ±S,, \A,\ ^ \Ay\ 

and corresponds to two massless Majorana fields with different velocities. This is a new, yet 
not considered possibility. 
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Conditions (|^) give four possible distributions of the exchange couphngs Ax,Ay,Bx and 
By on lattice bonds, each distribution breaking both U(l) spin rotational symmetry and 
translational invariance (below we indicate only the XF-part of the corresponding Hamil- 
tonians): 

H, = 2{A + B) Y.{SISI^, + SlSl^,) + 2{A - B) Y^i-^TiSlSl^, - SlSl^,), 

n n 

H, = 2iA - B) Y.{-ir{SnSl^i + SlSl^,) + 2{A + B) Y.{SISI_,, - SlSl^,), 

n n 

n n 
n n 

All these cases are related to each other by nonlocal unitary transformations. Transforming 
SI {-lySl, (-1)"^^, one finds that 

H2iA, B- J„ ^ HM. B; - J., -j',), H,{A, B; J„ j'^) ^ H,{A, B; - J„ -j'^). 

On the other hand, under — >• —5"^, 5"^ — > —5*^ only at n = 4j + 1 and n = 4j + 2, 

Hs{A, B- J„ J'^) ^ HM. B; - J„ -j',). 

In what follows, we shall concentrate on Hi. Choosing then Az = Bz, we arrive at the 
XXZ spin chain with XY bond-alternating (staggered) anisotropy, given by Eq.(0). It is 
equivalent to the following Hamiltonian describing two critical but different quantum Ising 
chains coupled by a self-dual interchain interaction: 

1 N/2 

^ i=i 

The model (^ has the following symmetry properties which will be used below: 
1) Under a global 7i/2 rotation SI -SI 

H{J,6,Jz)^H{J,-6,Jz). (8) 
7 



2) Under a staggered transformation, (— 1)"S'^, (— 1)""'"^S'^ 

H{J,5,J,)^H{-J,5,J,). (9) 

3) A nonlocal transformation to a 7r/2- twisted reference frame, 

S:^-S: (n = 4j + l,4j + 2), 
Sy^-Sy (n = 4j + 2,4j + 3), 

- (10) 

leads to the important symmetry under interchange of J and S: 

H{J,S,J,)^H{5,J,-J,). (11) 

III. XY CHAIN WITH BOND- ALTERNATING ANISOTROPY 

In this section we set Jz = and consider properties of the XY chain with a bond- 
alternating anisotropy: 

N 

H = Ho + Hs = T.[Jis:is^n+i + syAi) - 5i-ins:s:^, - sisy^,)]. (12) 

n=l 

Due to the symmetry property (P), we can choose without loss of generality J > 0. The 
Jordan- Wigner (JW) transformation 

n—l 2 
Sn = «n eXp(-Z7r a/fli). =■ ■= (^n(^n ' ' (13) 

i=l ^ 

represents Hamiltonian (|12D as a quadratic form of spinless fermion operators 

H= ^ ^[J(a+a„+i + a+_^ia„) - 5(-l)"(a+a++i + a„+ia„)], (14) 

describing a tight-binding half-filled band and a Cooper pairing of the fermions on neighbor- 
ing sites with a sign-alternating amplitude. The latter circumstance leads to commutativity 
of Hq and Hg which, in turn, indicates the absence of a gap in the spectrum at arbitrary 
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values of J and 6. This is in agreement with the equivalent representation of model in 
terms of two decoupled quantum Ising chains at criticality, as given by Eq.(0) at A = 0. 
Representation implies that the excitation spectrum consists of two branches which, 
in the continuum limit, correspond to two massless Majorana fields with group velocities 
v± = \J±6\ (we set the lattice constant a = 1). 

To understand the effect of the staggered XY anisotropy on the ground state properties 
of the model (p!4|), we rewrite the Hamiltonian in momentum representation 



H = Ho + Hs=Yl mata, + -A(fc)(a+atfc + h.c.)], (15) 

|fc|<7r ^ 

where 

e{k) = —e{k + tt) = Jcos /c, A(/c) = — A(A; + -n) = 5 cos k. 

We see that Hs describes Cooper pairing of JW fermions with total momentum equal tt. 
Mapping Hamiltonian (plSf) onto reduced Brillouin zone, |g| < 7r/2 



kl<V2 



and introducing a Nambu 2-spinor 

i^q = 



one rewrites H in the form 



H = Eo+ E < Kg) + A(g)ri] 7^,, (16) 

\q\<7r/2 

where r" are the Pauli matrices, and 

Eo = -T.<l) = -— (17) 

q 

is the ground state energy of Hq. 

Eq.([T6|) explains the origin of commutativity of Hq and Hs- The Hamiltonian of the 
"pure" XY model, Hq, with global U(l) spin rotation symmetry in real space, posesses a 
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local SU(2) pseudospin symmetry in momentum space0. Generators of this SU(2) group 
areEl 

'^," = ^^,+^>. {ci = x,y,z). (18) 

For a given g, A(g) then appears as a "magnetic" field oriented along the x-axis in the local 
pseudospin space. 

Diagonalization of H reduces to a 7r/2-rotation ipq — > exp[— (7r/2) J^] ipq, under which 

H = Eq + Hexcj 

H,x,= [E+{q)a^a, + E_{q)l3p,l (19) 

|9|<7r/2 

where new fermion operators 
and 

E±(g) = (J±5)cosg (21) 

At I (5 1 < J E±{q) > at all |g| < 7i/2. Therefore in Eq.(p!9D Hexc is positive definite 
and represents the excitation energy of the system, E±{q) thus being energies of the single- 
fermion excitations. The ground state \vac > is defined via relations aq\vac >= (3q\vac >= 
and coincides with that of Hq. The local pseudospin SU(2) symmetry is unbroken: 
J^\vac >= 0. 

Thus, at \6\ < J, the alternating XY anisotropy has no effect on the ground state of 
the system. This can be understood from the following simple arguments. Start with the 
ground state of the pure XY model Hq, 

\vac >= Y\_ '^fc |0 > 

7r/2<|fc|<7r 

and then consider Hg as a perturbation. Note that if a single-fermion state \k > is occupied, 
the state \7t — k > is necessarily empty, or vice versa. As a result, the Cooper pairing with 
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total momentum vr cannot occur in the ground state, if S is small enough; so the ground 
state is that of the XY model. This means that there is no off-diagonal long-range order 
(ODLRO) in the ground state 

< vac\a^a^_i^\vac > = 0, 
< vac\S:S:^, - SlS'^+,\vac > = 0, {\S\ < J), (22) 

implying that spin rotational U(l) symmetry and translational invariance remain intact. 
However, perturbation Hs reveals at the level of excitations. Consider one-particle and one- 
hole excitations \q >p= a^\vac >, jvr — g >h= aTT-q\vac > with |g| < tt/2. At 5 = these 
states are degenerate: Ep{q) = Efi{n — q) = Jcosg > (J > 0). The degeneracy is removed 
by Hs, but only at g 7^ ^kp, since the pairing amplitude A(g) has nodes at the two Fermi 
points. Therefore, a gap does not open; instead the doubly degenerate spectrum splits into 
two gapless branches with different Fermi velocities. 

The above picture persists at all \6\ < J. However, at \6\ > J one of the two excitation 
energies (|2l|) becomes negative, and the ground state changes. Since [Hq, Hs] = 0, this 
simply occurs due to level crossing: one of the former excited states becomes the ground 
state. As a result, the symmetry of the vacuum is changed. Choosing, e.g., 6 > J > 0, 
and making a particle- hole transformation (3q /3^, one finds that the ground state is 
now determined by aq\vac >= (3q\vac >= 0, the ground state energy and Hamiltonian of 
excitations being changed, respectively, by 

Eo = -EA(g) = -^, (23) 

i7e.c = E [E+{qH»i + E-iq)P^Pq] (24) 

|g|<7r 

with 

E±{q) = {5 ± J) COS q>0. (25) 

Notice that, at \S\ > J , Hq does not contribute to the ground state energy: 
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< vac\SlSl^, + SlSl^,\vac >= 0. (26) 

The local pseudospin SU(2) symmetry is now broken down to U(l), < vac\J^\vac >= —1/2, 
implying that the ground state is characterized by a nonzero ODLRO which breaks global 
U(l) and translational symmetries: 

1 flP 1 1 

° l^E(-l)" < v^clS^S:^, - SyS'^_,,\vac >= -9(6' - j') sgn 6, (27) 



N dS N^' ' I n n+i n n+i, ^ 

where 9{x) is the step function. Note that all results concerning the case |5| > J could have 
also been obtained by using the 5 J interchange symmetry, Eq. (^). 

Using the two-chain representation (ffl) of the spin operators, one finds that, at all values 



of J and (5, the local magnetization of the model (|T2D vanishes: < S*" >= 0. This simply 
follows from the fact that, since the a and r quantum Ising chains are decoupled and self- 
dual (critical), the Z2 x Z2 symmetry related to transformations (7| — — cr|, Tj — rj 
remains unbroken. Absence of local magnetization means that time reversal symmetry is 
always preserved. Preservation of time reversal symmetry together with breakdown of spin 
rotational U(l) symmetry allows to conclude that, at \6\ > J, the ground state represents a 
spin-nematic (SN) phasM, with the order parameter given by symmetric tensor (P?]). The 
SN state is doubly {Z2) degenerate; the two SN phases at 5 > J and 6 < — J transform to 
each other under transformation (|^). 

Let us now consider spin correlations in the ground state of the XY-like {\6\ < J) and 
SN {\6\ > J) phases. According to the two-chain representation of the Hamiltonian 
and multiplicative ar form (^) of the spin operators, the correlation functions factorize into 
products of independent contributions of each quantum Ising chain. At large space-time 
separations, the cr and (r) contributions are determined by gapless excitations, described in 
terms of (right- and left-moving) Majorana fermions with velocities f+ and f_, respectively. 
We shall first consider the case \6\ < J. Using the Majorana fermion representation of the 
quantum Ising chain0 and passing to a continuum limit, one can show that the longitudinal 
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{zz) spin correlation function is bilinear in the Majorana single-fermion Green's functions: 

1 / 1 / 1 



Att'^ \X + V^t X — V+t J \X + V-t X — v_t 
1 [1 - + [1 + 



(28) 



where x = n — m. This result can also be obtained by using the JW representation ([13|) for 
together with Bogoliubov transformation (pO]) to quasiparticle operators aq and The 
transverse spin correlation function is a product of two Ising correlation functions: 

< ^:W^^(0) > = < Syit)Sl{0) > - (29) 

At equal times, formulas (^) and (^) coincide with known results for the "pure" XY 
model0 (5 = 0). 

Now we consider the case \6\ > J. For simplicity, we shall restrict ourselves by equal 



time correlations. As follows from the 5 ^ J symmetry transformations ([T0|) , ([Tl]), in the 



SN phases the antiferromagnetic 2;2;-correlations (^) transform to the ferromagnetic ones: 

1 [1 - (-1)"-™] 

< ^n^m >vac= JT^ i f5 (30) 

/TT^ \n — my 

In the XY-like phase, the power-law decay (|2^) of the antiferromagnetic z2;-correlations 



leads to the well known singular response of the system to a staggered magnetic field along 
the z-axis, logarithmically divergent at T ^ 0. (In terms of the JW fermions, this prop- 
erty appears as a logarithmic charge-density-wave instability of a noninteracting ID Fermi 
system). Similarly, in the SN phases, ferromagnetic zz-correlations ( ppj ) give rise to a log- 
arithmically divergent response to a homogeneous magnetic field. The ferromagnetic and 
antiferromagnetic susceptibilities are finite at \5\ < J and \5\ > J, respectively. 
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At the boundaries 5 = ±J between the XY-hke and SN phases, the model shows both 
ferromagnetic and antiferromagnetic instabihties. 

Applying the 6 ^ J transformations (0) to (^), we find that, at \6\ > J , the trans- 
verse spin correlations show the same power-law decay as in the XY-like phase, but are 
characterized by a different type of local ordering: 

(-1)^ X X (-1)' 

< S2jS2j+2l >~ |2^|l/2' < ^2jS2j+2l+l >~ |2/ _^ ]^|l/2' ^^^"^ 

< SljSlj^2i >~ |2^|i/2' < 'S'|j'S'|j+2i+i >~ pipYpTi- (^^) 
Eqs.([3l|) and (|32D describe TTii and tilt periodic structures for the xx and yy spin correla- 



tions, respectively. The origin of these correlations is easily seen from (|T2D. At 5 > J the x- 
components of neighboring spins are coupled ferromagnetically on even bonds < 2j, 2j + 1 > 
and antiferromagnetically on odd bonds < 2j — l,2j >, whereas for the y-components the 
picture is opposite (at 5 < — J one has simply to interchange x- and y-components). 

The transverse spin correlations (|3TD, (|3^) locally describe a classical ground state of 
Hamiltonian Hs. Such a state represents a periodic structure, with a unit cell consisting of 
four spins = xcompm + ysiny^m, m = 1,2,3,4, where angles are given by ipi = ip, 
(p2 = n - (p, (ps = TT + (p, (p4 = -if at 6 > 0, and (pi = (p, (p2 = -(p, (ps = n + (p, 
ip4 = 7T — ip at 6 < 0, ip being an arbitrary angle. These configurations can be viewed as two 
interpenetrating antiferromagnetic sublattices whose staggered magnetizations transform to 
each other under refiection about the y-axis (5 > 0), or a;- axis (5 < 0), as required by the 
symmetry of Hs- Apparently, these states can be obtained from a classical Neel ground 
state of the XY model Hq by means of transformations ([TT|) and (H). Furthermore, the 
nonlocal transformation ( p!OD maps the continuous symmetry of Hq under uniform U(l) spin 
rotations onto the symmetry of Hs under continuous rotations of the two sublattices by 
the same angle but in opposite directions. As a result, the classical ground state of Hs is 
continuously degenerate under varying the angle ip. 
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In the quantum case of interest, the corresponding gapless "phason" mode destroys single- 
spin long-range order precisely in the same way as it occurs in the XY model. However, 
despite the loss of local magnetization, the survival of transverse correlations (0), ( |32D gives 
rize to SN ordering of the system, described by two-spin tensor (pT]). This is similar to the 
picture described by Chandra, Coleman and LarkinS for 2d frustrated antiferromagnets, 
where short-range correlations with local twist structure lead to the formation of a SN 
state. In the latter case, the order parameter is a parity-breaking antisymmetric tensor 
T(x, x') =< S(x) X S(x') >, corresponding to the p-type SN phasei. The difference between 
this and our cases is the absence of local twist structure in the ground state of our model at 
\5\ > J, (T(x, x') = 0). Instead we have the above described local spin structure generated 



by the staggered XY anisotropy, which results in a SN order parameter (pTj) with a symmetric 
tensor form. 

To conclude this section, let us briefly discuss the role of finite temperature. As it was 
already pointed out, the staggered XY anisotropy does not affect the ground state of the 
model at 6 < J, but reveals in splitting the excitation spectrum in two gapless branches 
with different Fermi velocities. One of interesting manifestations of this splitting is the 
appearance of the temperature-induced SN long-range order in the XY phase. Using the 
diagonalized form (|19D of the Hamiltonian, one finds that 

1 IE 

nx qx QV QV ___ -L 



, qx qx qy qy ^_ I 

Here E± are the average thermal energies of the a- and /5-quasiparticles. At low tempera- 
tures, T <^ J ±5, E±/N = (7r/24)(TV(^ ± 5)). Assuming for simplicity that 5 < J one 
finds 

< 5S5f,« >-< si_,si >^ [-2- + ^) + ^ (i) (f ^ 



< S^2,-iS2, >-< S^.S.^^i >- + 24 J2 ) I2\j)\j 
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The first term in the right-hand sides of these equations describes antiferromagnetic cor- 
relations in the XY model, slightly suppressed by thermal fluctuations. The second terms 
with opposite signs represent a combined effect of the temperature and XY bond-alternating 
anisotropy. These terms have the symmetry of the above described transverse spin corre- 
lations in the ground state of the SN phase and give rise to a finite value of the SN order 
parameter at \6\ < J 

< s:s':^, - syAi >= (-1)" (^) (33) 

The phenomenon we are dealing with here is an example of Villain's "order from 
disorder"!. Its manifestations in 2d frustrated magnets have been recently described in 
Refs.17,18. The 5-anisotropy cannot influence the ground state at 5 < J; however, it low- 
ers thermal excitation energy of the system by splitting the spectrum and repopulating the 
excited quasiparticles. This leads to the generation of new spin correlations, absent in the 
ground state, indicating the onset of SN order from thermal disorder. 

On increasing \6\, one of the two modes becomes very soft and nearly completely pop- 
ulated, when \ J — \6\ \ T <^ J + \6\. In this temperature range, the xx spin correlations 
on even (odd) bonds and yy correlations on odd (even) bonds get suppressed, and the am- 
plitude of the SN order parameter is close to l/27r. At \6\ J, the SN long-range order 
eventually penetrates into the ground state. 

IV. XXZ CHAIN WITH A WEAK BOND- ALTERNATING XY ANISOTROPY 

In this section, we extend the above considerations to the model (0) with a finite Jz and 
study effects of a weak bond-alternating XY anisotropy in the disordered gapless phase of 
the XXZ chain, assuming that \6\ ^ J, |A| < 1. The question we address here is: How 
off-diagonal perturbation Hg, tending to split velocities of fermionic excitations, reveals in a 
Luttinger spin hquidliiii, where orthogonality catastrophe suppresses single-fermion states 
in the low-energy part of the spectrum. Notice that the 6 ^ J symmetry (|TT]) allows to 
translate results of this section to the case of very strong 5-anisotropy, \6\ ^ J. 
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Being interested in the infrared properties of the model, we pass to a continuum descrip- 
tion developed for the XXZ spin chain by Luther and PeschelEl. Linearizing the spectrum 
of the JW fermions in the vicinity of two Fermi points ^kp = ±7r/2 and introducing two 
Fermi fields corresponding to right-moving [ipi) and left-moving (?/'2) particles 

an^(-0>i(a^)+^"^2(x), (34) 

one makes use of Abelian bosonizationlli'0 

: ?/'+2(a;)^i,2(a;) : = ^(9^v^i,2(a;) 

'71 



ipiA^) — exp[±iV^ V'i,2(a;)], (35) 

to describe the low-energy spin excitations in terms of a scalar field theory 

1 r 4A 
Hxxz = ^Vf / dx[. Ii\x) : +(1 + — ) : {dM^)f :] + Hu- (36) 

Here v^(x) and n(x) are the scalar field and its conjugate momentum, respectively. The 
term 

Hu--aJ dx cos{VlQ^<^{x)) (37) 



originates from Umklapp scattering of the JW fermiongH2llH3. At A > 1, these processes drive 
the system to a strong-coupling, massive Neel phase, but are irrelevant in the disordered 
phase, |A| < 1. Using continuum representation (^4|), the staggered XY anisotropy Hs 
transforms to a pairing term 

Hs = J dx[: iljt{x)^tix + a) : - : ^/'^(x)^/'^ (x + a) : -h.c] (38) 

Bosonizing (^) by means of ( |35D and performing canonical transformation of the field and 
momentum 

<^(x) ^ ^0(a;), U{x)^^P{x), PP = 4n, (39) 
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one arrives at the following self-dual field-theoretical model 

H = [ dx[^ (: P2(x) : + : (^^.^(x))^ :) + — cos/50(x) cos/30(x)]. (40) 

Here u is the renormalized velocity, a is a cutoff parameter, and is a field dual to 
defined as dx(f){x) = P{x). The correct parametrization of f3 is provided by the Bethe-ansatz 
solution of the XXZ modeS: 

co^.^A (41) 

Amazingly, the model ( ^OD is similar to the one recently proposed to describe two weakly 
coupled spinless Luttinger chainsi. The difference between the two models stems from 
different values of the conformal spin Sc = (3(3/271 of the perturbation: = 1 for the two- 
chain model, and Sc = 2 for the present one. As a result, the two-chain and present models 
belong to different universality classes. It was shown in Ref.8 that the self-dual theory 
(^) can be equivalently reformulated in terms of a 2d Coulomb gas of charge-monopole 
composites. The latter is given by the grand partition function: 

oo ^2n ^ 2n ^2^, / . \ 

Z/Z^ = E TTTj / n ^ E exp ^ a,a,k, +2KJ2 <^^^jhJ + 2z ^ (42) 

n=0 \'^^)- 3=\ " V i<3 i<j i^j I 

where is the partition function of the unperturbed, Gaussian part of the model. Eq.(^2p 
descibes a 2d system of classical particles with coordinates Xj = {uTj,Xj), (0 < Tj < 1/T), 
z = 6a /nu being the fugacity. Each particle carries an "electric" charge aj = ±1 and 
"magnetic" charge, or monopole, dj = ±1. The four- component system of charge-monopole 
composites is neutral with respect to each kind of charges. The charges and monopoles 
interact with logarithmic potential = /n(|xi — Xjl/a) amongst themselves, and there is 
also a statistical Aharonov-Bohm phase ipij = arctan((xj —Xj)/u{Ti — Tj)) which couples the 
charges and monopoles belonging to different composites. The parameters K = (3'^/Att and 
K = 1/K are naturally interpreted as inverse dimensionless temperatures for the charges 
and monopoles, respectively. 

The perturbation in ( |1D| ) has critical dimension D = K + K > 2. So, in the XY- 
like, gapless phase of the XXZ chain, the staggered XY anisotropy is marginal at the XY 
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point A = and irrelevant at A 7^ 0. The renormalized amplitude Seff{T), and hence the 
difference between the velocities of two Majorana fermionic excitations, scales down to zero 
according to the power law 6eff{T) ~ 5{aT/u)^~'^. This is analogous to the confinement 
regime in the two-chain modeli, in which Luttinger-liquid effects suppress single-particle 
interchain hopping in the infrared limilHi'i. 

At finite temperatures, thermally induced SN ordering will take place. However, in 
the Luttinger spin-liquid regime (A 7^ 0) , the SN order parameter increases with the 
temperature more slowly than in a noninteracting Fermi gas. Eg . (^3|) . Using Coulomb gas 
representation (|4^) and calculating the second-order correction to the free energy of the 
XXZ model, one finds that 

The nonuniversal power-law temperature dependence of the SN order parameter reflects the 
well-known infrared catastrophe in a ID Luttinger liquid that eliminates single- fermion states 
from the low-energy part of the spectrum. The r.h.s. of (^) should then be understood 
as ~ T^lNLLiT)]^, where Nil{u) ~ NQ{\u\a / u)^"""^ is the single-fermion density of states 
vanishing in the zero-energy limit. 

The suppression of the XY anisotropy in the lowest order in 6 does not really mean 
that Hg is a totally irrelevant pertubation. The usual criterium of relevance, based on the 
comparison of the critical dimension of a pertubation with space-time dimension 2, is not 
applicable here, since the pertubation has a nonzero conformal spin@. This is seen from the 
Coulomb gas representation (^21) . The suppression of the XY staggered anisotropy in the 
first order in 6 originates from pairing of the charge- monopole composites with zero total 
"electric" and "magnetic" charges. On the other hand, binding of composites in pairs with 
total "magnetic" or "electric" charge ±2 gives rise to two operators 

0{x) = cos2/50(x), 0(x) = cos2/30(x) (44) 

which are absent in the original model (|40|) , but are generated upon renormalization in the 
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effective Hamiltonian in the second order in (^i'^S. The operators (Q) have zero conformal 
spin; so their relevance can be estabhshed by the usual criterium. 



Clearly, 0{x) is the Umklapp operator (see (|37|)), relevant only at |A| > 1. 0(x) is a 
new operator whose critical dimension is 4K. It becomes marginal at the point = 1/2, or 
equivalently, A = Therefore, in the region — < A < 1, where both operators 0{x) 

and 0{x) are irrelevant, the model (|40|) preserves basic features of the gapless Luttinger- 
liquid phase of the XXZ chain, slightly modified by thermally induced SN ordering (^31) . This 
is a pure realization of Anderson's confinementi, i.e. stability of the Luttinger liquid against 
single-fermion off-diagonal perturbations, as opposed to the case of two coupled Luttinger 
chains, where the generation of two-particle interchain correlations drives inevitably the 
system away form the Luttinger-liquid fixed point at arbitrary in-chain interaction (see, e.g. 
Ref.24 and references therein). 

The point A = — is expected to be a Kosterlitz-Thouless transition point. The nature 
of the operator 0{x) and, accordingly, the type of ordering in a massive phase which occurs 
at A < — can be elucidated by means of the effective sine-Gordon model written in terms 
of the dual field 

Heff = I dx[j {P\x) + {d,~^{x)f) - ^cos2^0(x)], (45) 
together with the continuum representation of spin operators in terms of scalar fields (p and 

^ ^d.m + A(-l)"sin/?0(x), (46) 

SI ^ /x(-l)"cos^/30(a;), SI ^ /i(-l)" sin i/30(a;) (47) 

In the above equations, m is a positive parameter proportional to and A and are 

numerical constants. As follows from (^5]), the onset of a strong-coupling regime in the 
model ( |iO| ) results in ordering of the dual field, (f). 

Let us look at the symmetry properties of the model (9) under translations and spin 
rotations about the z-axisH. The spin operators in (10) and (11) are invariant under — »• 
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(/)+27r//5, (/) — > (/)+47r//3. This is the general symmetry of any spin modeL Under a translation 
by one lattice spacing (T^), (p + 7r//3, (p (f) + 2n/p. Under a U(l) rotation about the z- 
axis by an angle a (Ua), (p does not change, while 0+2a//3. Using these transformations 
of the fields and 0, one can check that the perturbation term in the original continuum 
model ( ^OD changes its sign under and Ut,/2, thus reflecting breakdown of translational 
symmetry (bond alternation) and the 7r/2 spin rotation symmetry (XY-anisotropy) . We 
then notice that, in the effective model P5|), each of these two symmetries is recovered; 
however the operator cos 2/30 changes its sign under 11-,,/^. We therefore conclude that 
the effective perturbation, generated in the second order in S, represents a translationally 
invariant, two-axis (Z4) anisotropy. Its lattice version is given by a four-spin term 

cos 2/30 ^ [S^S^j^i — 5'^5'^+l)(>S'^+2'5'j^+3 ~ Sn+2Sn+?i) 



fx \ 



( Qx cy _i_ Qy \( Qy 1 cy 

Wn'-'n+l "T "Jn^n+lJ Wn+2^n+3 "r '-'n+2^r! 

In the strong-coupling phase (—1 < A < —1/^/2) the Z4-symmetry is broken. The field 
acquires a nonzero vacuum expectation value, determined by one of the four degenerate 
minima = 0, 7r//5, 27r//?, Stt//? of the "potential" — mcos2/?0 (within the main periodicity 
interval Att/ (3). As follows from (0), this leads to a Neel ordering of the spins, either along 
the X-axis, 

(-ir<S:>~±M, {-IT<SI>=Q (48) 
at = 0, 27r//3, or along the y-axis, 

(-!)"< ^:>=0, {-lr<Sl>=^±^^ (49) 

at = 7r//3,37r//5. 

As we already mentioned, the results obtained for a weak bond-alternating anisotropy 
can be applied to the opposite limiting case \5\ ^ J, using the 6 ^ J symmetry (pH]), (p!0[). 



In the region —1 < A' < l/\/2, where A' = Jz/\5\, we find two gapless SN phases (for 
each sign of 6), with low-temperature properties of a Luttinger liquid. Applying nonlocal 



21 



transformations (0) to Eqs.(^8[) and (^9[), we find that at A' = l/\/2 a Kosterlitz-Thouless 
transition takes place to an ordered phase, characterized by the ttii periodic long-range 
alignment of the spins, either along the x- or y-direction in spin space. 

V. CONCLUSIONS 

We have shown in this paper that, in the S* = 1/2 quantum spin chains, the gapless, 
Luttinger spin-liquid state with a power-law decay of the correlation functions is not ex- 
hausted by systems possessing translational and spin rotation symmetries. Specifically, this 
type of low-energy behavior also characterizes the XXZ spin chain with the staggered XY 
anisotropy, exhibiting XY-like and spin-nematic phases with a gapless excitation spectrum. 
The model possesses a number of interesting properties. We have shown that massless Ma- 
jorana fermions with different velocities, being elementary excitations in the noninteracting 
case of the XY bond-alternating chain, fail to be observable (confinement) in the Luttinger- 
liquid regime, when interaction caused by a finite 22;— anisotropy A is switched on. Another 
interesting features of the model are order from disorder and Kosterlitz-Thouless transitions 
to massive ordered phases, driven by the 5— anisotropy. 

Although our description was restricted by limiting cases ^ J and \5\ J, it is clear 
that the gapless XY-like and SN phases occupy large domains in the parameter space of 
the model. However, to understand the phase diagram of the system in more detail and, 
in particular, describe a transition from the XY-like phase to the SN phases on changing 
5 with A kept finite, the region \5\ ~ J should be considered. This region is not accesible 
by bosonization method. As follows from the two-chain representation of the original 
model (^, in this region one has to consider an interacting system of "light" and "heavy" 
Majorana fermions with strongly different velocities (f+ ^ This problem resembles the 
Kondo-lattice one and deserves a special analysis. We hope to return to this and related 
questions in the future. 
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